Exact summation of vertex corrections to the penetration depth 

in d-wave superconductors 
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A variety of experiments suggest that in the cuprates, the low-energy superconducting quasi- 
particles undergo forward scattering from extended impurity potentials. We argue that when such 
potentials dominate the scattering, the penetration depth may be computed in a simple zero-angle 
scattering approximation (ZSA), in which the vertex corrections to the Meissner effect may be 
summed exactly. We find a remarkably simple relationship between the normal fluid density and 
the quasiparticle density of states of the disordered system which holds for every realization of the 
disorder. We expect this result to be relevant to the afe-plane penetration depth in high-purity single 
crystals of underdoped YBCO. 

PACS numbers: 74.25.Nf,74.72-h 



Measurements of the in-plane penetration depth A have 
been central to elucidating the nature of superconductiv- 
ity in the high-temperature copper-oxide superconduc- 
tors. Equally important has been the theoretical under- 
standing of impurity scattering in d-wave superconduc- 
tors required to interpret these measurements. The tem- 
perature dependence of A(T) provided the first strong 
evidencc"'^ for unconventional pairing symmetry. Subse- 
quently, the effect on A(T) of isotropic elastic scattering 
from point-like in-plane impurities in a superconductor 
with dx^-y^ pairing symmetry was investigatedi^ This 
picture of impurity scattering accounted for the observed 
crossover from a quadratic low-temperature behavior in 
A~^(T) to the T- linear behavior at higher temperatures 
characteristic of a gap function with nodes. 

More recently, it has been suggested in different ex- 
perimental contexts that quasiparticles near the Fermi 
surface undergo forward scattering. Forward scatter- 
ing from impurities has been invoked^ to account for a 
component of the single-particle scattering rate observed 
in angle-resolved photoemission (ARPES) spectra which 
apparently contributes negligibly to the normal state re- 
sistivity. Such scattering may arise from a smoothly 
varying in-plane potential due to the poorly screened"* 
Coulomb fields of disordered dopant oxygen ions residing 
between the Cu02 layers. This potential may account 
for the nanoscale electronic inhomogeneity observed in 
scanning-tunneling spectroscopy (STS) studies of opti- 
mally to ovcrdopcd Bi2Sr2CaCu2 08+a;i^ A flurry of re- 
cent activity^- '•^■^■iSi is aimed in part at synthesizing the 
evidence for forward scattering in this compound ob- 
tained from ARPES, STS, and thermal and electrical 
conductivities. 

At the same time, measurements of the microwave 
conductivity a{Lo,T) in YBCO present their own puz- 
zles. High quality YBa2Cu3 06.993 crystals have shown 
a Drude-like phenomenology in the conductivity at low 
frequencies and temperaturesjii In a nodal quasiparti- 
cle picture of low-energy charge transpor^i^^, this implies 
significant frequency-independent part of the scattering 
rate, inconsistent with scattering by point defects. It has 



been showni^ that extended linear defects such as twin 
boundary remnants can produce such behavior. How- 
ever, precise bolometric measurements^^ subsequently 
have shown deviations from Drude behavior at the lowest 
temperatures and frequencies as well as an unexpected 
Lo/{T + T{)) scaling. Violation of the Wiedemann- Franz 
law and the universal conductivity limit in cuprates 
have also fueled theoretical interest in quasiparticle self- 
energjii& and transportiii^ in the presence of extended 
impurities. 

Finally, extended disorder potentials have surfaced in 
a recent attempfeiS to explain in a unified theory the dop- 
ing and temperature dependence of both the c-axis and 
a6-plane penetration depth in underdoped YBCO. The 
theory models the incoherent interlayer hopping by scat- 
tering in-plane states by momenta of order A as they 
tunnel between the layers. Good agreement with recent 
c-axis data has been obtainedSP with h/K — 120A, or 
about 25 lattice spacings. The effect of such interlayer 
disorder on in-plane transport is similar to that of in- 
plane scattering potentials which extend over distances 
of order h/K. Presumably, such potentials arise from 
the aforementioned disorder of interlayer oxygen dopants, 
which, unlike BSCCO, may dominate the scattering rate 
due to the very low cation disorder in high-purity YBCO 
crystals. 

In this paper, we argue that the nodal structure of 
a d-wave superconductor presents an unusual situation 
in which even mildly extended potentials lead to strong 
forward scattering. The simplest approximation allow- 
ing for a sensible calculation of the superfluid density is 
the zero-angle scattering approximation (ZSA): the im- 
purity potential is unable to modify upon scattering the 
charge current carried by a quasiparticle. Extended po- 
tentials induce elastic scattering by small momenta, im- 
plying that, except via rare internodal scattering events, 
nodal quasiparticles scatter exclusively within the same 
node. However, the quasiparticle velocity relevant to 
charge current varies much more slowly in the nodal re- 
gion than does the energy, as it arises from the hare band- 
structure. The ZSA is thus expected to be valid when 
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extended potentials dominate the scattering rate. This 
is likely the situation in high-purity YBCO single crys- 
tals, in which smooth potentials due to interlayer dopants 
dominate the scattering. 

The main result presented here is that under such cir- 
cumstances the normal fluid density nn is connected to 
the single-particle density of states per site per spin M{E) 
by the simple relation 

[ dEN{E)[-nE)i (1) 

m* \ a / J 

where vp is the magnitude of the quasiparticle velocity at 
the node and / is the Fermi-Dirac distribution. Since the 
electron density is the sum of the normal- and superfluid 
densities, this quantity dictates the temperature depen- 
dence of the penetration depth. It is remarkable that 
([T|). which is a trivial result for the clean superconduc- 
tor, continues to hold in the presence of extended disorder 
potentials, albeit with M{E) renormalized by the disor- 
der. Since it is valid for every realization of the disorder, 
it circumvents the considerable complications of vertex 
corrections^ which arise upon disorder-averaging two- 
particle correlation functions. Furthermore, ^ does not 
require any particular model of the disorder (e.g. Born 



limit, ladder approximation,) though it reveals nothing 
about U{E) itself. Thus if the ZSA is valid, iP) is a 
powerful consequence allowing a model-independent in- 
terpretation of the normal fluid density. 

The limitation of this approximation is that, as 
in the clean superconductor, the dissipative part of 
conductivity^ (Ji(uj) becomes a (5-function at zero fre- 
quency. This is to be expected since zero-angle scatter- 
ing cannot degrade the charge current. However, as in a 
clean superconductor at T 7^ 0, a fraction of the super- 
fluid is converted into normal fluid, which may be reliably 
computed despite the fact that the distribution of ai [uj) 
is unphysical. In the ZSA, this normal fluid behaves as 
a perfect metal, having infinite d.c. conductivity but 
making no contribution to the Meissner effect. In real- 
ity, small-angle (and rare internodal) scattering degrades 
this perfect metal and results in a dissipative conductiv- 
ity peak of finite width. However, the ZSA can provide 
no information about the frequency distribution of ai (w). 
We assume that such scattering is sufficiently weak that it 
primarily affects only the frequency distribution of ui [uj) 
and not the integral, so that Eq. remains valid. 

We consider the general Hamiltonian H — Hk + Hgc: 
with 



r, CJ 



h.c. , 



i +Cr+ya,tCr, I.) 



(2) 



h.c. 



It describes, at the mean-field level, a c?-wave supercon- 
ductor on square-lattice layers with in-plane lattice con- 
stant a and interlayer spacing d. The interlayer hopping 
(which may model the incoherent c-axis transport as 
in Ref. and the on-site potential Vj- are assumed to 
be weak relative to t and have a zero average value. The 
order parameter Ar is defined on lattice bonds and may 
also vary weakly (compared to its average value) due to 
the disorder potentials. These quantities are assumed to 
vary slowly along the in-plane directions over the length 
scale h/ K but may vary randomly from layer to layer. 

It will be convenient to work with the Nambu spinor 
^I'r = (^r, p j,)"^- The Hamiltonian (0) may be diag- 
onalized by a Bogliubov transformation to H = Ea + 
I]n,<Te"7i,o-7n,<T where £^0 is the ground state energy 
and e„ labels the (positive) quasiparticle excitation en- 
ergies. The 7 operators obey the usual fermion com- 
mutation relations up to a renormalization constant: 

{7r cT' 7r' ct'} ~ N6r^r'5(j,a' whcrc N in a finite system 
is the number of lattice sites. The transformation can 
be compactly represented using a two-component wave- 



function u'^"''(r) — (M„(r), w„(r))-^. Introducing the anti- 
symmetric tensor e, we write (with the sums over tensor 
indices impficit hereafter,) 7^ | = Sr(^r)"^""H'') ^'^'^ 

7ri,i = I]r£aa'(*J)a'uL"''*(r). The wavefunctions Sat- 
isfy the orthogonality relations 

l^ui")(r)ui"')*(r)=<5^„S 

1 < (3) 

^^e.,u(r)(r)uf )(r) = 0. 

r 

With this notation, the in-plane components of cur- 
rent are J,^y){v) = it{^U^(^^^j! ^ " *J*r+x(y)a)- 1^ 
the disordered system, nn/m* — T\xx{Q) where the 
long- wavelength polarization function is Tlf^u{i^) = 

^ Er, r' So dr e»-- {%J^{r, r) Mr', 0)). It is also useful 
to introduce a vertex function which describes the cou- 
pling of long-wavelength radiation to electrons at lattice 
sites ri and r2, 

Aa;(y)(l"l,l"2) = ii((5ri,r2+x(y)a - (^rs, ri+x(y)a)- (4) 
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Importantly, the Hamiltonian Q describes a non- 
interacting system, and we may apply Wick's theorem 
to evaluate U. We assume that (J) = in the unper- 
turbed system, in which case 11 may be expressed 



T ^Tr [G(r2 ,r[;i:y + iuj)G{r'^,ri;ii^)] 



(5) 



in terms of the matrix Green function G(r, r';ij/) = 
-/(f dr e""^ <r^^'r(r)^'|:,(0)>. We will make use of 
the spectral function A{r,r';uj) = (27ri)~^ [G'^(r, r'; oj) — 
G^(r,r';aj)] which has the representation 



71 

+ 5{uj + e„)eaa/e/3;5/u|^"^(r')u^7^ *(r) 



(6) 



We now state the ZSA precisely, establish its validity, 
and derive (Q). We define 



(7) 



whose trace and integral over w gives 2 Ax- The quantity 
Qx may be regarded as a decomposition of A^ in terms 
of quasiparticle energies oj. Eq. (0 for the polarization 
becomes 



N ^ 



dE dE 



, f{E) - f{E') 
lu - [E' - E) (8) 



Tr[Q,(ri,r2;£;)Q,(r2,ri;i;')]. 

The ZSA is implemented in the disordered system by 
approximating Qx by a form which is valid for low-lying 
quasiparticle energies E, E' . 

To find such an approximation, it is instructive to view 
Aa;(ri, r2) as the kernel of an integral operator, which we 
call the "vertex operator" . To illustrate, we rewrite Q 
as 



Ax(ri,r2) = 



(2^ 



2isin(fcxa)e''^-""i (e*'"''^)*. (9) 



The vertex operator thus has eigenvalues 2tsm{kxa) and 
plane-wave eigenvectors. On the other hand, the quasi- 
particle wavefunctions of the clean system (i.e. — 0, 
= 0, and Ar constant) are also plane waves. The 
(5-function behavior of conductivity in the clean sys- 
tem arises precisely because the vertex function and the 
Hamiltonian are diagonal in the same basis, i.e. the ver- 
tex operator commutes with H . In the disordered sys- 
tem, H is diagonalized by the wavefunctions u^"^ and 
thus does not commute with the vertex operator, mak- 
ing cTi {lj) regular. 

A low-energy effective vertex operator for the clean 
system is obtained by restricting k in to the vicinity 



of a node q. Under this restriction, the eigenvalue is 
weakly k-dependent. Thus we make the approximation 



2tsin(A;j.a) -^vp j a\pi^ 



(10) 



the sign depending on the node. The leading k- 
dependent term (fc^ — qx)/m*a leads in the clean sys- 
tem to a correction to whose relative size we esti- 
mate as (2'r]T / m* Vp)^ where 77 > 1 is the anisotropy 
of the linearized nodal dispersion. For the particular 
bandstructure of (0), rn*Vp = 4i sing^^atangajO, which is 
nominally of order At but diverges rapidly as half-filling 
is approached. Furthermore, r]Tc <SC 4t in underdoped 
cuprates, making an excellent approximation. 

The approximate vertex operator arising from Hl()|l is, 
within the subspace of a single nodal region, simply a 
constant multiple of the identity operator. This powerful 
simplification may be applied to the disordered system 
given two assumptions. We first assume that the low- 
energy wavefunctions in the disordered system are linear 
combinations of low-energy wavefunctions in the clean 
system. We thus neglect the contribution of low-energy 
states which may arise from the mixing of anti-nodal 
stateSiiS. Since the anti-nodes are connected by wavevec- 
tors of order a~^, such low-energy states should consti- 
tute a negligible portion of the spectral density when 
Aa <C 1. Second, we assume that the extended impuri- 
ties produce negligible internodal scattering. Degenerate 
plane-wave states from different nodes will hybridize and 
split by an energy of order \V'\, where V' is a typical 
matrix element for scattering between nodes. Assuming 
the Fourier spectrum of Vr to be exponential with width 
A, we estimate \V'\ ~ Vbe"'*''!/^ where Sq is the rela- 
tive momentum between nodes and Vq is the RMS value 
of the potential Already \V'\ < Vq when A^^ is 

comparable to just a few lattice spacings. We assume 
that relevant experiments are in the regime T ^ \V\, so 
that these hybridized states are essentially degeneratei2& 
We may then choose a basis for the low-energy sector of 
the disordered system in which each state resides near a 
single node. 

Thus, in the regime \V'\ T ^ 4</77, it becomes 
trivial to express the vertex operator in the quasiparticle 
basis of the disordered system. We denote the eigenvalues 
as Ai"^ = ±VF/aV2, the sign depending on the node of 
state n. Then Q becomes 



<5(c.-e„)uf)(ri)u(")*(r2) (11) 
d{uj + e„) ■Hr2)u(7)*(ri)l, 



which may be regarded as the precise statement of the 
ZSA. Since the time reversal transformation on wavefunc- 
tions is ui"''(r) — > Eaa'U^"' *(r), the second term in ifTT)) 
corresponds to the time-reversal of the first and acquires 
a sign change. The orthogonality relations allow us 
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to evaluate 

^ Tr[Q,(ri,r2;i?)g,(r2,ri;ii;')] = (12) 
i (^) ' - E')^Y.^S{E e„) + 5{E + 

n 

Substituting into ||SJ) leads directly to main result J^l, 
with the density of states per site per spin M given by 

n 

We note that (^1 can also be obtained^i using the di- 
agram methods^ for disorder-averaged correlation func- 
tions. In this context, zero-angle scattering at the nodes 
allows vertex correction diagrams to be summed exactly 
and expressed in terms of the single-particle Green func- 
tion. This simplification is formally analogous to the di- 
agrammatic Ward-Takahashi identity-?? of quantum elec- 
trodynamics, by which the longitudinal component of the 
interacting vertex function may be expressed in terms of 
the interacting Green function. In our context, however, 
has nothing to do with electromagnetic gauge invari- 
ance, but arises instead from the fact that i) the Hamil- 
tonian Q is noninteracting and ii) the effective vertex 
function for the low-energy states has trivial structure, 
as we have argued above. 

In ordinary metals forward scattering leads to a re- 
duction of the transport scattering rate 1 /rtr compared 
to the single-article scattering rate l/r. The result pre- 
sented above can be viewed as an extreme example of 



this effect in a system of nodal fermions where l/r is 
finite but l/rtr is reduced essentially to zero. This illus- 
trates the singular nature of the vertex corrections that 
have been summed to obtain It is also interesting 
to note that the superfluid density, ordinarily considered 
a transport quantity, depends only on the single-particle 
scattering rate in this case. 



We have derived the remarkably simple result 
which relates the normal fluid density to the quasiparti- 
cle density of states in a system disordered with extended 
impurities. The normal fluid density may be computed in 
the ZSA in spite of the fact that, as in a clean supercon- 
ductor, the real part of conductivity becomes a 5-function 
at zero frequency. This derivation proceeds directly from 
the quasiparticle wavefunctions u'"'* and energies e„ of 
Hamiltonian ((2Jl, showing that holds for individual 
realizations of the disorder and thus circumvents the con- 
siderable complications of vertex corrections. It therefore 
allows for a model-independent inversion of n„(T) from 
experiment to obtain the density of states. To theoreti- 
cally predict N{E) and cri(tt'), however, requires a more 
specific, model-dependent calculation. Eq. |Q may be 
applicable to high-purity single crystals of underdoped 
YBCO, in which slowly varying potentials from interlayer 
dopant disorder dominate the low energy scattering. 
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Hereafter we use a-\{LL}) and "conductivity" to denote only 
the normal fluid component. 

Generally, A and t^ will have spatial variation as well as 
V . In this case, Vb characterizes the largest of the three 
fluctuations. 

Beyond the ZSA, the splitting of these states governs the 
broadening of the (5-function in conductivity. 



